Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



X-ray Line Emission from Supernova Ejecta Fragments 

A. M. Bykov 



(N 
O 
O 



in 



> 
in 

(N 
(N 

in 
o 

O 

Oh 
6 



A.F. Ioffe Institute for Physics and Technology, St. Petersburg, 194021, Russia 
Received 4 March 2002 / Accepted 26 April 2002 

Abstract. We show that fast moving isolated fragments of a supernova ejecta composed of heavy elements should 
be sources of K a X-ray line emission of the SN nuclear-processed products. Supersonic motion of the knots in the 
intercloud medium will result in a bow-shock/knot-shock structure creation. Fast nonthermal particles accelerated 
by Fermi mechanism in the MHD collisionless shocks diffuse through a cold metallic knot, suffering from Coulomb 
losses and producing the X-ray emission. We modeled the X-ray emission from a fast moving knot of a mass 
M ~ 10 -3 Mq , containing about ~ 10 -4 Mq of any metal impurities like Si, S, Ar, Ca, Fe. The accelerated 
electron distribution was simulated using the kinetic description of charged particles interacting with strong MHD 
shock. We accounted for nonlinear effects of shock modification by the nonthermal particles pressure. The K a line 
emission is most prominent for the knots propagating through dense molecular clouds. The bow shock should be 
a radiative wave with a prominent infrared and optical emission. In that case the X-ray line spectrum of an ejecta 
fragment is dominated by the low ionization states of the ions with the metal line luminosities reaching L x > 
10 31 erg s _1 . High resolution XMM-Newton and Chandra observations are able to detect the K Q and Kp line 
emission from the knots at distances up to a few kpcs. The knots propagating through tenuous interstellar matter 
are of much fainter surface brightness but long-lived. The line spectra with higher ionization states of the ions are 
expected in that case. Compact dense knots could be opaque for some X-ray lines and that is important for their 
abundances interpretation. The ensemble of unresolved knots of galactic supernovae can contribute substantially 
to the iron line emission observed from the Galactic Center region and the Galactic ridge. 

Key words. Acceleration of particles; Radiation mechanisms: non-thermal; ISM: clouds; ISM: individual object: 
IC 443; ISM: supernova remnants; X-rays: diffuse bacground 



^ 1. Introduction 



Physical mechanisms of energy release and nucleosynthesis 
products ejection in supernova (SN) events are of funda- 
mental importance for astrophysics. Optical and UV stud- 
ies of the structure of SN remnants (SNRs) have revealed 
a complex metal ejecta structure with the presence of iso- 
lated high- velocity fragments of SN ejecta interacting with 
surrounding media. The most prominent manifestations 
of this phenomena are the fast moving knots (FMKs) ob- 
served in some young "oxygen-rich" SNRs in the Galaxy: 
Cas A (e.g. Chevalier & Kirshner 1979; Fesen et al. 2002), 
Puppis A (Winkler & Kirshner 1985), G292.0+1.8, in the 
LMC N132D and IE 0102.2-7219 in the SMC. 

Optical FMKs in Cas A are showing very high abun- 
dances in O-burning and Si-group elements with high 
Doppler velocities. They have a broad velocity distribu- 
tion around 6,000 km s _1 , while bright radio knots are 
slower (e.g. Bell 1977). The optically emitting mass of the 
observed knots is only about lO~ 4 A/ (Raymond 1984). 

Ejecta-dominated X-ray emission from Cas A was ob- 
served with Chandra (Hughes et al. 2000; Hwang et al. 
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2001) and XMM-Newton (Bleeker et al. 2001; Willingale 
et al. 2002). Willingale et al. (2002) found an emis- 
sion component that is probably due to heating of ejecta 
clumps sweeping up the ambient gas. They also note that 
the observed Fe-K line emission is confined to two large 
clumps. Recent Chandra images of G292.0+1.8 (Park et 
al. 2001) revealed a complex structure with multiple knots 
and filaments on angular scales down to the instrumental 
resolution. 

X-ray knots were discovered with ROSAT in the 
Vela SNR by Aschenbach et al. (1995) as 'shrapnels', 
boomerang structures outside of the main shell. High- 
resolution Chandra observations of shrapnel A reported 
by Miyata et al. 2001 revealed a head-tail structure of 
the apparent size 8'. 4 x 4'.1 (0.6 x 0.3 pc at 250 pc dis- 
tance). They estimated the gas pressure in the head to be 
roughly ten times higher than that in the tail and T e ~ 
0.5 keV. The oxygen abundance was 0.34+q;q| of the solar 
value while that of Si was times of solar. 

A localized region of 6.4 keV emission indicating the 
presence of Fe XVII or lower ionization states was found in 
the supernova remnant RCW 86 with ASCA instrument 
by Vink et al. 1997. 
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High resolution Chandra observation of the Galactic 
Central region reported by Wang et al. (2002) revealed 
that the observed He-like Fe K a emission is largely due 
to the discrete X-ray source population. Apart from the 
large-scale thin plasma emission in the Galactic Center re- 
gion Bamba et al. (2002) have also observed with Chandra 
many small clumps of emission lines from neutral (6.4 
keV) to helium- like (6.7 keV), with intermediate line en- 
ergies between 6.5-6.7 keV. The problem of the origin of 
the observed large scale X-ray emission from the Galactic 
ridge requires a careful study of possible classes of abun- 
dant hard X-ray sources with L x ~ 10 29 erg s _1 (e.g. 
Tanaka et al. 1999). In this respect modeling of an en- 
semble of X-ray point sources associated with fast moving 
SN ejecta fragments seems to be necessary. 

In the paper we present a model of X-ray line produc- 
tion in a supersonically moving ejecta fragment. The line 
emission is due to if-shell ionization by nonthermal par- 
ticles accelerated by the bow shock and then propagating 
through a cold metal-rich clump. The optical depth effects 
are important in the model and they are accounted for. 
Modern arcsec resolution instruments such as that aboard 
XMM-Newton and Chandra are sensitive to photons up to 
10 keV with a good spectral resolution providing unique 
possibilities to study X-ray line emission of the ejecta frag- 
ments. These observations would allow to study stellar 
nucleosynthesis. 

2. Model 

Consider a fast-moving metal-rich ejecta fragment. 
Supersonic motion of the fragment (of velocity v k , ra- 
dius R and mass M) in the ambient medium will result in 
a bow-shock/internal-shock structure creation (Chevalier 
1975; Sutherland & Dopita 1995). The bow shock has the 
velocity > v k and it is propagating in an ambient medium 
of the standard cosmic composition. The bow shock cre- 
ates a hot high pressure gas and also accelerates nonther- 
mal particles. The high pressure gas in the head of the 
fast moving ejecta fragment could drive an internal shock 
of velocity w; s into the fragment if w; s exceeds the sound 
speed in the cold matter: 

V is « V k (pa/pk) 1 / 2 OC V k p^R^M-^ (1) 

Here p a and pk are the ambient gas and the dense frag- 
ment densities, respectively. We shall also use the density 
contrast \ = Pk/Pa- The characteristic size R of 10 16 - 
10 17 cm is consistent with that resolved by optical ob- 
servations of Cas A and Puppis A while large clumps 
were inferred by Chandra observations of Vela shrapnel 
A by Miyata et al. 2001. For an ejecta fragment of radius 
i?i6 ~ 3 (where Rie is the radius measured in 10 16 cm) 
and mass M ~ 10 _3 M Q the density ratio p a /pk = X -1 is 
about 4-10 -4 for an ambient medium of number density 
~ 0.1 cm~ 3 . However the ratio could be as high as > 0.1 
for the ejecta fragments moving through a dense ambient 
matter of a molecular cloud. 



There is a number of factors limiting the fast ejecta 
fragment lifetime. A fast moving knot is decelerating due 
to the interaction with the ambient gas. The deceleration 
time of a knot can be estimated as Mv\j p & v\S ~ 

xR/v-k, where S is the effective crossection of the bow 
shock. However, the fast knot hydrodynamical crushing 
effects could have much shorter timescale, providing r c ~ 
X 1/2 i?/«k ~ R/v is oc M 1 /2(i? . p a )- 1 /2 u -i ( e .g. Wang 
& Chevalier 2001 and references therein). That implies 
Tc/Td ~ x~ 1/2 oc i? 3 / 2 p a /2 M" 1 / 2 . The knot ablation pro- 
cess with surface erosion could also be a limiting factor. 

The hydrodynamical estimation of the inner shock ve- 
locity Vi s given by the Eq.(l) assumes an efficient conver- 
sion of the bow shock ram pressure to the knot internal 
shock. That is a good approximation if the effects of en- 
ergy loss due to nonthermal particle acceleration and ra- 
diative cooling of high pressure gas behind a bow shock 
are negligible. The latter is the case unless the knot is in- 
teracting with dense matter (like molecular cloud clumps). 

The Alfvcnic Mach number of both the bow shock and 
the inner shock (if Vi s is determined by the Eq.(l)) is 

M & = v k (47rp a ) 1 / 2 /B w 460 v s n x J 2 / B_ 6 , (2) 

where n a is the ionized ambient gas number density mea- 
sured in cm~ 3 , _B_6 is the local magnetic field just before 
the shock measured in pG and v$ is the flow velocity in 
10 8 cm s _1 . Note here that the Alfvenic Mach number 
of the inner shock depends on the ambient density but on 
the magnetic field inside the knot. The sonic Mach num- 
ber for a shock propagating in a plasma of the standard 
abundance is 

M s ^8Hv 8 -pva + Zei)]- 1 / 2 , (3) 

where T4 is the plasma ion temperature measured in 10 
K and / e i = T e /Tj. In general, for the processes in the pre- 
cursor and viscous velocity jump of an MHD collisionless 
shock wave no equilibration between electrons and ions 
should be assumed (e.g. Raymond, 2001). Both cases of 
/ e i <C 1 and / e i > 1 could be relevant to the system under 
consideration. The photoionized gas at the radiative pre- 
cursor of a fast shock is assumed to have T ~ (1 — 2) • 10 4 
K and / ei > 1 (e.g. Shull & McKee 1979). On the other 
hand, plasma heating due to MHD waves dissipation in 
the vicinity of the viscous subshock heats mostly the ions, 
providing / ei < 1. 

2.1. The effect of nonthermal particles on a strong 
bow shock 

For strong collisionless shocks in a magnetized plasma 
the nonthermal particle acceleration effect is expected to 
be efficient and a significant fraction of ram pressure is 
transferred to high energy particles (ions for nonrelativis- 
tic shocks). The shock transition of a strong shock of the 
total Mach number Mt 3> 1 is broadened because of the 
upstream gas deceleration by nonthermal particle pres- 
sure gradient ahead of the viscous gas subshock of a mod- 
est Mach number .Msub < 3 (see for a review Blandford 
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& Eichler 1987; Jones & Ellison 1991; Malkov & Drury 
2001). 

The total compression ratio R t of a strong MHD shock 
modified by an efficient nonthermal particle acceleration 
depends on the energy flux carried out by escaping non- 
thermal particles and the effective adiabatic exponent (e.g. 
Malkov & Drury 2001). On the other hand the distribu- 
tion function of nonthermal particles and the bulk flow 
profile in the shock upstream region are sensitive to the 
total compression ratio R t . Thus, the exact calculation 
of the escape flux can be performed only in fully nonlin- 
ear kinetic simulations. Nevertheless, some approximate 
iterative approach can be used to make the distribution 
function consistent with the shock compression. The sub- 
shock is the standard gas viscous shock with the com- 
pression ratio: R s ^ h = ( 7g + l)M s 2 ub / [( 7g - l)M 2 suh + 2]. 
The downstream ion temperature behind the modi- 
fied shock structure can be calculated if Rt and i? S ub are 
known for the shock of the given velocity and it could 
be very different from that obtained with the standard 
single-fluid Rankine-Hugoniot law. The electron tempera- 
ture just downstream of the shock is expected to be lower 
than T^l 

An exact modeling of the collisionless shock struc- 
ture with nonthermal particle acceleration effect requires 
a nonperturbative selfconsistent description of a multi- 
component and multi-scale system including strong MHD 
turbulence dynamics. That modeling is not feasible at the 
moment. We shall use some simplified description of a 
strong shock with some apropriate parameterization of 
governing parameters which are: (a) the nonthermal par- 
ticle diffusion coefficients; (b) the ion injection rate; (c) 
the maximum momentum of accelerated particles. The 
gas heating mechanism due to MHD waves dissipation in 
the shock precursor must also be specified in that simpli- 
fied approach, providing the connection between the total 
compression and the gas subshock Mach number A^sub- 
Assuming that the main heating mechanism of the gas in 
the precursor region is due to Alfven waves dissipation, 
Berczhko & Ellison (1999) obtained a simple relation for 
the total compression of the shock R t ~ l.BM^ 8 , while 
the subshock had A4 sub ~ 3. The relation is valid un- 
der the condition of A4 2 3> A4 a in the far upstream flow. 
They also noted that the subshock compression ratio is 
not too sensitive to the ion injection rate if the rate ex- 
ceeds 5 • 10~ 5 to provide the shock modification by the 
accelerated ions. We obtain then the ion temperature just 
behind the modified strong shock (measured in 10 6 K) 

T 6 (2) « 0.32 • 4>(M suh ) ■ v 5 k i 4 n- 3 ^B% 4 (l + f ci )-\ (4) 

where #M Bub ) = [2 7g A4 2 ub -( 7g -l)]/[( 7g -l)A4 2 ub + 2], 
and (f)(3) ps 3.7 for 7g = 5/3. Then the electron and ion 
temperatures equilibrate in the postshock region. The ini- 
tial electron temperature in the shock downstream de- 
pends also on the collisionless electron heating (see e.g. 
Bykov & Uvarov, 1999; Laming 2001). The turbulent gas 
heating due to acoustic instability wave dissipation could 



be more efficient however, if the acoustic instability devel- 
ops (Malkov & Drury 2001). 

The maximal momenta of the accelerated ions are lim- 
ited by MHD wave damping due to ion-neutral collisions 
in the case of a fragment interacting with dense molecu- 
lar cloud. The effect was accounted for following the work 
by Drury et al. (1996). The maximal momenta of the 
ions are < 10 2 GeV/c in that case. In the case of a frag- 
ment propagating through tenuous medium the maximal 
momentum is limited by the finite scale of strong scatter- 
ing region with typical values of the maximal momentum 
> 10 2 GeV/c. 

It should be noted that we still have no direct observa- 
tional evidences for the cosmic ray ions accelerated in SNR 
shocks. At the same time there are direct evidences for 
ion and electron acceleration by strong cosmic-ray mod- 
ified shocks in the interplanetary medium (e.g. Terasawa 
et al. 1999). In the interplanetary shock of velocity ~ 375 
km s _1 (in the solar wind frame) studied by Terasawa et 
al. (1999) with the GEO TAIL satellite the partial pres- 
sure of accelerated nonthermal particles was found to be 
in a rough balance with the thermal proton pressure in the 
shock downstream. That implies a high efficiency of non- 
thermal ion production. The diffusion coefficients for both 
electrons and ions were estimated to be <~ 10 18 cm 2 s _1 in 
the shock vicinity that is at least two order of magnitude 
smaller than the typical interplanetary values (Shimada 
et al. 1999). For SNR shocks Dorfi (2000) estimated the 
ion acceleration efficiency to be 0.24 for some reasonable 
parameters. 

The observed synchrotron emission from SN shells is 
considered for a long time as a direct evidence for relativis- 
tic electron acceleration. The efficiency of transformation 
of the kinetic energy of a nonrelativistic shock bulk flow 
into the nonthermal electron population is not expected to 
be very high, being below a few percent. However even a 
percent range efficiency of the nonthermal electron accel- 
eration would provide an observable effect of K Q photon 
emission from metallic knots irradiated by high energy 
electrons. Bulk heating and ionization of the knot interior 
due to deep penetration of energetic ions and electrons 
could also be important for enhanced ablation. 

2.2. Nonthermal electron acceleration 

To simulate the spectra of nonthermal electrons in the 
fragment body we solve a kinetic equation for the distri- 
bution function of the electrons, accounting for injection, 
diffusive transport, advection, first and second order Fermi 
acceleration and synchrotron and the Coulomb losses. 
Electron kinetics in supercritical collisionless shocks was 
modeled by Bykov & Uvarov (1999). They showed that 
strong MHD fluctuations generated by kinetic instabilities 
of ions are responsible for heating and pre-acceleration of 
nonthermal electrons on a very fine scale, of the order of 
several hundreds of inertial lengths, in the vicinity of the 
viscous jump (subshock) of a collisionless shock. We have 
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described the details of the code, diffusion model and loss 
functions used elsewhere (Bykov & Uvarov 1999; Bykov 
et al. 2000). 

There are three zones in the one-dimensional model: 
the pre-shock region (I), the shock transition region (II), 
and the post-shock flow (III), from where nonthermal 
emission from shock-accelerated particles originates. In or- 
der to calculate the spectra of nonthermal electrons in 
these regions, we used the following kinetic equation for 
the nearly- isotropic distribution function Ni (z, p) (i = I - 
III): 

, , , d 2 N l (z,p) 1 d 2r ^ . , dN t 

Hp) — ^ + _ _ ADi{p) _ + Li{p) Ni] _ 

This Fokkcr-Planck-type equation takes into account 
diffusion and advection [bulk velocity Ui{z)\ of electrons 
in phase space due to interactions with MHD waves and 
the large-scale MHD flow. Here z is the coordinate along 
the shock normal. Li(p) is the momentum loss rate of an 
electron due to Coulomb collisions in a partially ionized 
plasma and synchrotron/inverse Compton radiation. The 
momentum diffusion coefficient D(p) is responsible for the 
second order Fermi acceleration, and fcj (p) is the fast par- 
ticle spatial diffusion coefficient. For low-energy electrons 
the Coulomb and ionization losses are important every- 
where except for the narrow shock transition region (II), 
where acceleration is fast enough to overcome the losses 
and where nonthermal electron injection occurs. 

To accelerate the electrons injected in the shock transi- 
tion region to relativistic energies MHD turbulence should 
fill the acceleration region. As we know from the inter- 
planetary shock observations, turbulence generation is a 
generic property of collisionless MHD shocks. We used the 
parameterization of the momentum dependent diffusion 
coefficient ki (p) of an electron of momentum p, described 
in details by Bykov et al. (2000). The diffusion coefficients 
are momentum dependent with ki(p) oc k ■ vp above 1 
keV. The normalization value ko is the electron diffusion 
coefficient at 1 keV. 

The electron temperature in the (photo)ionized shock 
precursor was fixed to T c = 2 • 10 4 K for the fragments 
propagating in the dense clouds of the gas of number den- 
sity 7$> 1 cm~ 3 . The fragment velocity should exceed 
~ 250 km s _1 to provide enough UV radiation to com- 
plete H and He (to He + ) ionization in the precursor of 
the bow shock, if it is modified by the nonthermal par- 
ticle pressure. The velocity is higher than the value of 
~ 110 km s" 1 derived by Shull & Mckee (1979), because 
of the weaker gas heating by the modified MHD shock 
structure with a smooth precursor and a gas subshock of 
■Msub ~ 3. It should be noted however that even at lower 
shock velocities where the upstream photoionization is in- 
efficient, the preshock ionization could be provided by the 
fast particles accelerated at the collisionless shock. For the 
fragments propagating in hot tenuous plasma the electron 



temperature was assumed to be equal to the plasma tem- 
perature (i.e. / e i = 1). 

The energy spectrum of the nonthermal electrons in 
the bow shock of a metal-rich ejecta fragment is shaped 
by the joint action of the first and second order Fermi ac- 
celeration in a turbulent plasma with substantial Coulomb 
losses. The electron spectra arc calculated for different val- 
ues of the bow shock velocity assuming that the diffusion 
coefficient is parameterized as in Bykov et al. (2000) with 
different values of fco for 1 keV electron for different en- 
vironments. The minimal value of the diffusion coefficient 
of an electron is k^ in « 6.6 • 10 16 B~^(E /keV) cm 2 s" 1 . 
Note that the minimal value is only achievable if strong 
fluctuations of magnetic field are present at the electron 
gyro-radii scale. For slightly superthermal electrons the 
fluctuations cannot be generated by the ions. Whistler 
type fluctuations could provide the required electron scale 
fluctuations (e.g. Levinson 1996). However in the case of 
a strong MHD turbulence the electron transport by the 
vortex type fluctuations of ion gyro-radii scale would pro- 
vide the effective diffusion coefficient much larger than 
k™ m (Bykov & Uvarov, 1999). We consider here rather a 
conservative case of that large diffusion coefficient. 

The electron acceleration model described above is 
one-dimensional. The transverse gas flow of velocity u± 
behind the bow shock front would affect the nonthermal 
particle distribution. The exact 2-D and 3-D modeling of 
the postshock flow with account of shock modification ef- 
fect is not feasible at the moment. However, the trans- 
verse flow effect on high energy electrons is only marginal 
if the particle diffusion time across the transverse flow 
is shorter than the transverse advection time scale i.e. 
6±/R ■ u±S±/k(p) < 1. Here S± is the width of the trans- 
verse gas flow behind the bow shock (typically it is a rel- 
atively small fraction of the fragment radius R). In our 
case the condition is satisfied for 100 keV - MeV regime 
electrons dominating the K a lines production in the ejecta 
fragment body. We account for the presence of a postshock 
cooling layer of the shocked ambient matter situated in 
front of the metal rich fragment. The curves 1-3 in Fig. 1 
illustrate the effect of Coulomb losses on low energy non- 
thermal electrons diffusing through the postshock layer. 

2.3. Ambient conditions and the fragment state 

As a generic example we consider a fast moving fragment 
of ~ 1,000 km s _1 , mass M = 10~ 3 M Q and radius 
Ri6 = 3 dominated by oxygen and containing about <~ 
10~ 4 Mq of any metal impurities like Si, S, Ar, Ca, Fe. 
That fragment we refer as "the standard knot". 

The electron distribution was simulated here for two 
distinct representative cases of the ambient environment. 
The first one is an SN fragment propagating through dense 
ambient medium of number density ~ 10 3 cm~ 3 , relevant 
to the SN explosion in the vicinity of a molecular cloud. If 
the ambient magnetic field follows the scaling B a oc n 1 / 2 
(e.g. Hollenbach & McKee 1989) one may expect the mag- 



A. M. Bykov: X-ray Line Emission from Supernova Ejecta Fragments 



5 




Fig. 1. The nonthermal electron distribution function 
(N(p, x) — N(p,x)p 2 pt) as a function of dimcnsionless 
particle momentum p/pr (pt is the momentum of a ther- 
mal electron in the upstream region) . We show the distri- 
butions at different distances from the bow shock plane 
for different ambient media: (a) a fragment of velocity 
t>k = 3, 200 km s _1 in a tenuous medium; (b) an SN frag- 
ment of velocity = 1, 080 km s _1 in a molecular cloud 
(see Section 2.3 for details). The depths from the bow 
shock plane are 10 6 , 10 12 , 10 15 , 3- 10 16 cm for the curves 
1-4, respectively. 



netic field values of > 100 [iG in the dense ambient 
medium (c.f. Blitz ct al. 1993). The second case is a SN 
fragment propagating through tenuous medium of num- 
ber density ~ 0.1 cm~ 3 , temperature T > 10 4 K and 
magnetic field of /j,G range. Correspondingly, the diffusion 
coefficient normalization factors for keV regime electrons 
were k ~ 10 18 cm 2 s _1 in the dense medium case and 
ko Si 10 20 cm 2 s _1 for the tenuous medium. The gas ion- 
ization state in the bow shock downstream was simulated 
here in the way like that described in the paper of Bykov 
et al. (2000). 

The gas ionization state inside an oxygen rich fragment 
body was estimated following the model by Borkowski & 



Shull (1990). They studied the structure of a steady state 
radiative shock (inner shock) of velocities V[ s from 100 
to 170 km s _1 (c.f. Eq.(l)) in a pure oxygen gas with 
electron thermal conduction. The postshock temperature 
behind the 140 km s^ 1 shock with thermal conduction 
was found to be 270,000 K in the immediately postshock 
region and dropped to roughly 4,700 K after the oxygen 
depth about 2- 10 15 cm -2 in the downstream photoioniza- 
tion zone. The dominant ionization state of oxygen in the 
shock radiative precursor and in the immediate postshock 
region was Oil I while it was Oil in the downstream pho- 
toionization zone. Nonthermal electrons could also be an 
important factor to model the dominant ionization state 
of the ions in the fragment body (Porquet et al. 2001). 

In Fig. 1 the nonthermal parts of the electron distri- 
bution functions N(p,x) = N(p,x)p 2 pT are shown at the 
different depths of an SN fragment. Hcrepx is the momen- 
tum of a thermal electron in the far upstream region. Note 
that the upstream gas electron temperatures arc different 
for the two ambient media we displayed in the Fig. 1. In 
Fig. 1 we showed only nonthermal parts of the electron 
distributions because the electron temperatures in both 
presented cases are below those required for efficient K a 
lines production in the ejecta fragments. The threshold for 
K-shell ionization of Si is 1.8 keV (i.e. p/pt > 33 in Fig. 
1) while that for Fe is about 7.1 keV (i.e. p/pt > 64 in 
Fig. 1). Fast electrons accelerated in the MHD collision- 
less shock diffuse through the postshock layer and cold 
metallic knot suffering the Coulomb losses as it is clear 
from the Fig.l. They are producing the K a lines in the 
fragment body due to radiative transitions following the 
removal of the Is atomic electrons. 

2.4. K-shell ionization and X-ray line emissivity 

To calculate the production rate of characteristic X-ray 
line emission from the metal rich fragment irradiated by 
an intense flux of energetic electrons accelerated by an 
MHD shock wave the atomic inner shell ionization cross 
sections and the fluorescent yields are required. The ex- 
perimental electron-impact i-T-shell ionization crossections 
and the fluorescent yields were compiled recently by Liu 
et al. (2000). The standard Bethe theory of inner shell ion- 
ization (e.g. Powell 1976) was used to fit the data for non- 
rclativistic energies of ionizing electrons, from the thresh- 
old up to 50 keV and the relativistic formulae by Scofield 
(1978) above 50 keV. 

We fit the K-shell ionization cross section following the 
parameterization given by Scofield (1978) 

a l = ^-(b + b 1 + b 2 /e + b 3 /p 2 + b 4 ■ b/p + b 5 /p 4 ) (6) 

Here the crossection Oi is given as a function of the in- 
cident electron momentum p (measured in m e c), b = 
ln(p 2 ) — f3 2 , (3 — v/c. The fitting parameters A, b\ — 6 5 
were calculated by Scofield (1978) for the ions of Z > 
18. We corrected the fitting parameters to account for the 
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present laboratory measurements compiled by Liu et al. 
(2000). 

The logarithmic asymptotic of the crossection at rela- 
tivistic energies, given by the Eq.(6) is important for our 
modeling. One can see in the Fig.l (curves 3 and 4) that 
the electron spectra have rather flat maxima at mildly rel- 
ativistic energies (p/pr <^ 400) in the most of the volume 
of a dense metal- rich fragment. Thus, the K-shell ioniza- 
tions by relativistic electrons are substantial. 

The radiative decay of K-shell vacancy induced by fast 
electron results in K a (2p — * Is) and Kp (3p — > Is) lines 
production. The I(Kp)/ T(K a ) intensity ratio is about 
0.14 (e.g. Scofield 1974) for Fel-FelX and it is decreasing 
for higher iron ionization stages (e.g. Jakobs & Rozsnyai 
1986). The ratio is an increasing function of the atomic 
number Z. Note however, that the optical depths for Kg 
and K a lines could be different, especially for highly ion- 
ized iron, and that would affect the observed line ratio. 
Thus, observing the ratio I(Kp)/ I(K a ) one may constrain 
the source optical depth. Both I{K a ) and I(Kp) Fe lines 
were detected recently with Chandra from Sgr B2 giant 
molecular cloud by Murakami et al. (2001). The observed 
ratio is somewhat below 0.2 indicating some possible op- 
tical depth effect, but it is still marginally consistent with 
the value 0.14 expected for a transparent system. 

It is important to note here that shock-accelerated en- 
ergetic nuclei can also provide efficient K-shell ionization. 
The emission spectra from the decay of inner-shell ioniza- 
tions produced by a collision with an MeV-regime ion have 
multiple satellites (peaks) at energies higher than that of 
the electron-induced K a line (e.g. Garcia et al. 1973). We 
shall discuss the effect of ions in detail elsewhere. 

2.5. The optical depth effect 

The maximal column density of a metal (of atomic 
weight A) in a spherical ejecta fragment NmeL = 5.6 • 
10 22 (M^/Mq ) A^ 1 cmT 2 can be high enough 

to provide optical depths r > 1. We use below N^ A > = 

(l/2)7\4iax as the resonant scattering column density in 
the model assuming the shock situated at the fragment 
center. The resonant line scattering effect (see, for the re- 
view, Mewe 1990) can be important for the metal rich SN 
fragments. The optical depth at the line center of an ion 
due to the resonant line scattering was given by Kaastra 
& Mewe (1995) and can be expressed as: 

r a 1.34 / abs iV 1 ( 5 A) J B pli 1 ( J 4/T) 1 / 2 (l + 0.5^/T)- 1 / 2 . (7) 

where £" p h (in keV) is the line photon energy, / a b s is the 
absorption oscillator strength ,T is the ion temperature 
(in eV), iV-fjj is the ion column density (in 10 15 cm -2 ), 
wq is the gas micro-turbulence velocity (in 10 km s _1 ). 

To describe the radiative transfer for the line radiation 
produced inside the absorbing medium we used a simpli- 
fied escape probability approach. The mean escape prob- 
ability pf(r) can be approximated as Pf(r) m [1 + ar] _1 , 
where a « 0.43 for r < 50 (Kaastra & Mewe 1995). The 



Table 1. K-shell line luminosities of the fragment inter- 
acting with a molecular cloud 



Line(j 


1,080 


( km s 
1,620 


- 1 ) 
2,700 


'max 


O (0.54 keV) 


40.4 


104 


1638 


33,880 


Si (1.7 keV) 


1.3 


10 


48 


592 s 


Ar (2.9 keV) 


0.5 


4 


20 


272 fc 


Fe (6.4 keV) 


0.4 


3 


15 


78 6 


T m [ 10 7 K ] 


0.3 


0.5 


1.0 





a The luminosities are in 10 38 ph s" 1 
h The absorption depths can be applied only for the 
ionization states Si VI, Ar X, Fe XVIII and higher. 

photon absorption due to the resonant scattering of K- 
shell line is most important for the ions with an incom- 
plete L-shell (e.g. OI, Si VI, Ar X, Fe XVIII) and relatively 
low fluorescent yields. The resonant absorption is not ef- 
fective for the neutral or the low ionization stages of Fe, 
Ca, Ar, Si with completed L-shells where we applied the 
photoabsorption by Morrison & McCammon (1983). 

3. X-ray emission from an SN fragment 
interacting with a molecular cloud 

Core-collapsed SNe from massive progenitors are expected 
to be correlated with massive molecular clouds. The most 
spectacular manifestation of the phenomena are the SNRs 
in starburst galaxies like NGC 253, Arp 220, M82 etc. 
(see e.g. Chevalier & Fransson 2001). There are also star- 
forming molecular clouds Sgr A, Sgr B2 in the Galactic 
Center region with a strong CS emission from dense molec- 
ular gas (e.g. Blitz et al. 1993). To model such a case we 
considered a SN fragment of radius 3T0 16 cm and of oxy- 
gen mass of 10 -3 Mq , containing also ~ 10~ 4 M© of an 
impurity (Fe,Ar,Si). The fragment is propagating through 
a molecular clump of a number density 10 3 cm~ 3 with a 
velocity > 1,000 km s~ x . The magnetic field value in the 
cloud is about 100 pG and fco « 10 18 cm 2 s _1 at 1 keV. 

We calculated the local emissivities of K-shell lines 
through the fragment depth as well as the integrated line 
luminosities. Table 1 contains the luminosities of K Q lines 
of O, Si, Ar, Fe. These luminosities are not corrected for 
the optical depth effect. We have also given the maximal 
depths r max for resonant scattering calculated under the 
assumption that the column density of a given charge state 
of the ion A is equal to Nimh.. To obtain a real estimation 
one should correct the depth for the actual column den- 
sities of the ions with incomplete L-shells. Note that in 
Tables 1 and 2 we presented only the integral luminosi- 
ties of K Q complex. A prediction for the line shape de- 
pends on the details of the ion charge state profile inside 
the fragment. The exact modeling of the ionization profile 
is beyond the scope of the present paper. Nevertheless, 
one could see that in the case when the bow shock is 
radiative (see below) the dominant iron charge state in 



A. M. Bykov: X-ray Line Emission from Supernova Ejecta Fragments 



7 



the fragment becomes lower than Fe XVIII, providing the 
low resonant absorption depth and the line centroid to be 
close to 6.4 keV. We assumed in the Table 1 that the gas 
micro-turbulence velocity Wq — 1, and the ion tempera- 
ture T < 10 4 K in the fragment body. The ion tempera- 
tures T^ 2 ) presented in the Tables 1 and 2 were calculated 
for the position just after the viscous subshock. The ob- 
servable electron temperatures in the postshock relaxation 
region are somewhat lower. 

As seen in Table 1, an increase in the fragment ve- 
locity results in a strong enhancement of the line lumi- 
nosities. For the lower knot velocities the line emission 
drops down drastically because the Coulomb losses are 
dominating in that case. Note however that if the value of 
fco was less than 10 18 cm 2 s" 1 (at 1 keV) even the frag- 
ments of lower velocity (i>k < 1, 000 km s _1 ) could provide 
a substantial line luminosity. A fast fragment interacting 
with CS-emitting gas of density > 10 4 cm~ 3 could reach 
even higher line luminosity if fco ~ 10 17 cm 2 s _1 (at 1 
keV). Since the larger fragments (for a given knot mass) 
are more transparent for the K-shell lines of oxygen (and 
other elements) a fast SN fragment of R J> 3 • 10 17 cm, 
propagating through the inter-clump medium of density 
> 10 cm -3 , might have a prominent oxygen K-line. 

A fast fragment of a larger scale > 10 17 cm, of the 
same mass 10 -3 M© , entering the molecular clump will 
have p a / pk > 1 and would drive a strong shock into the 
metal-rich fragment. Such a fragment should be a source 
of gamma-ray lines and also light (and other spallogenic 
origin) elements produced by accelerated ion interactions 
with the metal-rich knot. It would appear as a bright tran- 
sient source. 

4. X-ray emission from an SN fragment in a low 
density medium 

Many of the isolated SN fragments are propagating 
through a low density hot environment. That concerns 
both a relatively low velocity fragments moving inside the 
forward shock radius of SN and fast velocity fragments of 
an SN exploding in a low-density environment of the num- 
ber density n xs 0.1 cm~ 3 , gas temperature T rs 2 • 10 4 K 
and the magnetic field value ~ 3 fiG. We first considered 
the same SN fragment as that in the dense medium, but 
in a wider velocity range 1, 000 < «k < 7, 000 km s _1 be- 
cause even very high velocity fragments are long-lived in 
the low density environment. The diffusion coefficient nor- 
malization at 1 keV was fixed to be fco ps 3 • 10 19 cm 2 s _1 . 
We summarize the simulated X-ray line luminosities (mea- 
sured in 10 36 photon s~ 4 ) in Table 2. One can see that the 
X-ray line luminosities L x < 10 29 erg s _1 (per 10 -4 M Q 
of Si, S, Ar, Ca, Fe) are predicted from an individual SN 
fragment of the scale 3T0 16 cm in a tenuous medium. In 
Table 2 we assumed that the gas micro-turbulence veloc- 
ity Wq = 1, and the ion temperature T < 10 4 K in the 
fragment body. 

However, simulations of somewhat larger > 10 17 cm 
fragments show that the luminosity corrected for the ab- 



Table 2. K-shell line from the fragment interacting with 
low-density gas 



Linejj 


1,600 


( km s 
3,200 


- 1 ) 
6,400 


Tmax 


O (0.5-0.6 keV) 


38.0 


66.5 


99.8 


33,880 


Si (1.7-1.8 keV) 


2.6 


4.5 


6.7 


592 5 


Ar (2.9-3.1 keV) 


1.9 


3.4 


5.0 


272 fc 


Fe (6.4-6.9 keV) 


0.8 


1.4 


2.1 


78" 


T (2) [10 7 K] 


1.2 


2.8 


6.6 





a The luminosities are in 10 36 ph s" 1 
h The absorption depths can be applied only for the 
ionization states Si VI, Ar X, Fe XVIII and higher. 



sorption is increasing to L x > 10 30 erg s _1 and even 
higher due to decreasing of the optical depth and Coulomb 
losses. It is important that the large scale fragments are 
much thinner (and hotter) providing a substantial amount 
of ions in high ionization states. These faint transparent 
fragments would contribute substantially to the observed 
diffuse X-ray line emission of highly ionized matter. Note 
that the fragment deceleration time is oc Mn~ l R~ 2 . That 
implies that the fragments in the old remnants could only 
be observed if they spent most of the time in the tenuous 
medium. 

4.1. A model of the shrapnel A in the Vela SNR 

We simulated the line emission from a fragment of a 
scale ~ 10 18 cm to model the Vela shrapnel A discov- 
ered by Aschenbach et al. (1995) and recently stud- 
ied with Chandra by Miyata et al. (2001). The oxygen- 
dominated fragment of mass M ~ 10~ 2 M Q and velocity 
v k ^ 10 8 cm s~ 4 would have the deceleration time about 
10,000 years in an ambient medium of n w 0.1 cm~ 3 . 
We found that the temperature behind the fragment bow 
shock dominated by nonthermal particles is about 0.5 keV. 
The silicon line at 1.8 keV would have the luminosity L x ~ 
10 30 erg s -1 if ~ 1O _3 M of Si is contained in the frag- 
ment and the oxygen line at 0.6 keV - L x > 10 31 erg s _1 . 
The resonant absorption depth of the Si K-shell line is 
Tmax ~ 0.5 while that of oxygen is r max ~ 30 assuming 
the gas micro-turbulence velocity wq — 10, and the ion 
temperature T > 10 6 K in the large diluted fragment. The 
mean escape probabilities Pf(r) are about 0.8 for Si and 
0.07 for oxygen. The optical depth effect could account for 
the apparent Si overabundance observed by Miyata et al. 
(2001). 

5. Discussion 

There are observational evidences for the clumpy structure 
of SN ejecta in the remnants of different types and ages. 
Observations are available for the core-collapsed SNRs 
(e.g. Cas A) and also for the remnants are thought to 
be of Type la (see e.g. the recent XMM-Newton observa- 
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tions of Tycho by Decourchelle et al. (2001)). Some of the 
SNRs are young (e.g. Cas A) and others are older (e.g. the 
Vela SNR). 

The observation of the ejecta clumps composition 
would provide a valuable test to study the details of SN 
phenomena. The simulations of SN nucleosynthesis (e.g. 
Woosley & Weaver 1995; Thielemann et al. , 1996) pre- 
dict the production of ~ O.1M of 56 Ni and 28 Si and 
> O.OO3M of 36 Ar in a core-collapsed SN. The metal 
distribution through the ejecta at the postcxplosion stage 
is still poorly known. There are conclusive evidences for 
large-scale macroscopical mixing of Ni in the SN 1987A 
ejecta (e.g. McCray 1993). However the modeling of the 
56 Fe distribution through the ejecta at the microscopic 
level is rather a complicated task. This is because of the 
energy deposition effect from 56 Ni decay (see e.g. Wang & 
Chevalier 2001; Blondin, Borkowski & Reynolds 2001 for 
recent discussions). The X-ray line emission in our model 
is less sensitive to the details of the element distribution in 
a fragment because of the high penetrating ability of en- 
ergetic electrons responsible for the line excitation. Even 
if the Fe (or other element) atoms are locked in the dust 
grains we would still have X-ray line emission excited by 
fast electrons. 

We show that the X-ray line luminosities and spec- 
tra are different for the SN ejecta fragments in a dense 
molecular cloud and for that in a tenuous medium. The 
efficiency of X-ray line production is higher in the frag- 
ments of lower ionization state. Strong MHD shocks could 
transform a sizable fraction of the kinetic energy into non- 
thermal particles thus reducing the heating of the shocked 
gas and increasing its compression. The column density 
of the shocked gas could reach ~ 10 21 cm~ 2 and the 
ion temperatures just behind the shock transition region 
is below 5-10 6 K. Thus, the fast fragments (of velocity 
v k is 1,500 km s _1 ) moving through molecular clouds of 
density > 10 3 cm~ 3 would drive the radiative bow shock 
wave. Hollenbach & McKee (1989) have modeled the spec- 
tra of the radiative shocks. They found rich spectra of H2 
emission as well as atomic fine-structure lines with strong 
[OI] (63 //to) [OIII], [Fe II], [C II] (158 //to), [Si II] (35 
Hm), [Ne II] (12.8 //to) lines. That implies a correlation 
of IR and optical emission of the fragment with the X-ray 
line emission dominated by the lines of relatively low ion 
charge states. The soft thermal X-ray continuum emission 
from the postshock gas of T < 10 7 K is not expected to be 
observable in that case. On the other hand the nonthermal 
brcmsstrahlung emission with the hard spectrum of a typ- 
ical photon indexes < 1.5 is predicted in our model. The 
continuum bremsstrahlung emission is rather sensitive to 
the ionization state of the fragment. Our simulations show 
that the X-ray line luminosity is < 10~ 4 of the kinetic 
energy dissipation rate in the fragment's bow shock. 

The lifetime of the X-ray line emitting fragments in 
a dense molecular cloud is expected to be a few times 
less than the fragment deceleration time. Such a fragment 
should appear as a nonthermal source showing a variabil- 
ity on a timescale of years. The variability is most im- 



portant for radio emission. Synchrotron radio emission at 
a level about 100 mJy arcsec -2 at 100 MHz could be ex- 
pected from the fragments interacting with a dense molec- 
ular cloud if the ambient magnetic field is about 0.1 mG. 
Radio emission from a fragment propagating through a 
tenuous medium is much fainter. However, MHD type in- 
stabilities may greatly enhance the magnetic fields along 
the fragment boundary providing substantial radio emis- 
sion, see Anderson et al. (1994); Jones et al. (1996). 

The X-ray line luminosity L x ~ 10 31 erg s~ 4 (per 
10~ 4 M Q of Si, S, Ar, Ca, Fe) is predicted from an in- 
dividual SN fragment in a molecular cloud. It could be 
detected from a few kpc distance with the current de- 
tectors aboard Chandra and XMM-Newton . An obvious 
candidate to study is IC 443: the SNR interacting with 
a molecular cloud. The recent XMM-Newton observation 
of IC 443 by Bocchino & Bykov (2001) has resolved the 
hard X-ray source 1SAX J0618. 0+2227 correlated with a 
molecular cloud into two sources one of which is generally 
consistent with the ejecta fragment interpretation. A ded- 
icated observation of 1SAX J0618.0+2227 is desirable to 
check the interpretation. 

An ensemble of unresolved SN fragments could con- 
tribute substantially to the observed diffuse iron line emis- 
sion. Diffuse iron line emission has been found by ASCA 
from the Galactic Centre region on a scale about a degree 
along the galactic plane (Koyama et al. 1996) and from a 
more extended region of the Galactic Ridge (e.g. Tanaka 
2002). He-like and H-like Fc-K lines analyzed by Tanaka 
(2002) are significantly broadened corresponding to a ve- 
locity dispersion of a few thousand km s~ 4 . The high 
dispersion can be naturally explained if the lines are due 
to unresolved faint SN fragments. Alternatively, Tanaka 
(2002) has suggested that the high dispersion is due to 
charge-exchange processes of low-energy cosmic rays pro- 
viding the high velocities of the line emitting ions. 

A molecular cloud irradiated by a hard X-ray source is 
expected to be observable in fluorescent X-ray lines as X- 
ray reflection nebulae - a new class of X-ray sources (e.g. 
Sunyaev et al. 1993; Sunyaev & Churazov 1998). The ob- 
served iron line emission from the Sgr B2 (Murakami et al. 
2001) and Sgr C (Murakami et al. 2001a) complexes was 
considered by the authors as the X-ray reflection nebulae 
irradiated by a bright source in the GC which was active 
in the past. Clumps of X-ray line emission of neutral and 
highly ionized iron of equivalent widths ~ keV and absorp- 
tion corrected 2-10 keV fluxes ~ 10~ 12 erg s _1 cm~ 2 was 
found by Chandra in the Galactic center region (Bamba 
et al. 2002). The photon indexes of 2-10 keV continuum 
emission - while not too constrained - are broadly consis- 
tent with the photon indexes < 1.5 predicted by the model 
of SN fragments emission considered above. If a substan- 
tial fraction of the iron ejected by an SN in a molecular 
cloud would reside in the fast moving fragments we could 
expect the iron line luminosity of L x ~ 10 34 erg s _1 per 
SN which is consistent with that observed from the Sgr C 
complex. The line emission from the Sgr B2 is somewhat 
brighter and would require contributions of more than one 
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SNc. Infrared line observations from radiative shocks that 
accompany SN fragments in molecular clouds would help 
to distinguish the SN fragments contributions. A correla- 
tion between SiO (J=2^0) emission morphology and 6.4 
keV Fe line emission on the Galactic Center large scale and 
also within Sgr A and B regions found recently by Martin- 
Pintado et al. (2000) could be consistent with the multiple 
SN fragments contribution. ISO observations of 18 molec- 
ular clouds in the Galactic Center region were performed 
by Martin-Pintado et al. (2000a). Towards most of the 
clouds they have detected the "ionized bubbles" with the 
fine structure line emission of ionized species: SIII, Nell, 
and in some cases also Nelll, Nil, NIII, OIII. The "bub- 
bles" could be relevant to the fragment-cloud interactions. 

Nonthermal hard X-ray continuum emission from the 
Galactic ridge requires a population of accelerated elec- 
trons (e.g. Valinia et al. 2000, Dogiel et al. 2002). Hard 
continuum of the SN fragments could contribute to the 
observed emission. The galactic diffuse emission mapping 
with forthcoming INTEGRAL mission will address the is- 
sue. 
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